home *** CD-ROM | disk | FTP | other *** search
/ Chip 2007 January, February, March & April / Chip-Cover-CD-2007-02.iso / Pakiet multimedia / Muzyka / Edytory sampli (probek dzwieku) / ZynAddSubFX_2.2.0 / Setup_ZynAddSubFX-2.2.0.exe / source code / DSP / FFTwrapper.C < prev    next >
C/C++ Source or Header  |  2005-03-14  |  3KB  |  100 lines
  1. /*
  2.   ZynAddSubFX - a software synthesizer
  3.  
  4.   FFTwrapper.c  -  A wrapper for Fast Fourier Transforms
  5.   Copyright (C) 2002-2005 Nasca Octavian Paul
  6.   Author: Nasca Octavian Paul
  7.  
  8.   This program is free software; you can redistribute it and/or modify
  9.   it under the terms of version 2 of the GNU General Public License 
  10.   as published by the Free Software Foundation.
  11.  
  12.   This program is distributed in the hope that it will be useful,
  13.   but WITHOUT ANY WARRANTY; without even the implied warranty of
  14.   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15.   GNU General Public License (version 2) for more details.
  16.  
  17.   You should have received a copy of the GNU General Public License (version 2)
  18.   along with this program; if not, write to the Free Software Foundation,
  19.   Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  20.  
  21. */
  22.  
  23. #include <math.h>
  24. #include "FFTwrapper.h"
  25.  
  26. FFTwrapper::FFTwrapper(int fftsize_){
  27.     fftsize=fftsize_;
  28.     tmpfftdata1=new fftw_real[fftsize];
  29.     tmpfftdata2=new fftw_real[fftsize];
  30. #ifdef FFTW_VERSION_2
  31.     planfftw=rfftw_create_plan(fftsize,FFTW_REAL_TO_COMPLEX,FFTW_ESTIMATE|FFTW_IN_PLACE);
  32.     planfftw_inv=rfftw_create_plan(fftsize,FFTW_COMPLEX_TO_REAL,FFTW_ESTIMATE|FFTW_IN_PLACE);
  33. #else
  34.     planfftw=fftw_plan_r2r_1d(fftsize,tmpfftdata1,tmpfftdata1,FFTW_R2HC,FFTW_ESTIMATE);
  35.     planfftw_inv=fftw_plan_r2r_1d(fftsize,tmpfftdata2,tmpfftdata2,FFTW_HC2R,FFTW_ESTIMATE);
  36. #endif
  37. };
  38.  
  39. FFTwrapper::~FFTwrapper(){
  40. #ifdef FFTW_VERSION_2
  41.     rfftw_destroy_plan(planfftw);
  42.     rfftw_destroy_plan(planfftw_inv);
  43. #else 
  44.     fftw_destroy_plan(planfftw);
  45.     fftw_destroy_plan(planfftw_inv);
  46. #endif
  47.  
  48.     delete (tmpfftdata1);
  49.     delete (tmpfftdata2);
  50. };
  51.  
  52. /*
  53.  * do the Fast Fourier Transform
  54.  */
  55. void FFTwrapper::smps2freqs(REALTYPE *smps,FFTFREQS freqs){
  56. #ifdef FFTW_VERSION_2
  57.     for (int i=0;i<fftsize;i++) tmpfftdata1[i]=smps[i];
  58.     rfftw_one(planfftw,tmpfftdata1,tmpfftdata2);
  59.     for (int i=0;i<fftsize/2;i++) {
  60.     freqs.c[i]=tmpfftdata2[i];
  61.     if (i!=0) freqs.s[i]=tmpfftdata2[fftsize-i];
  62.     };
  63. #else
  64.     for (int i=0;i<fftsize;i++) tmpfftdata1[i]=smps[i];
  65.     fftw_execute(planfftw);
  66.     for (int i=0;i<fftsize/2;i++) {
  67.     freqs.c[i]=tmpfftdata1[i];
  68.     if (i!=0) freqs.s[i]=tmpfftdata1[fftsize-i];
  69.     };
  70. #endif
  71.     tmpfftdata2[fftsize/2]=0.0;
  72. };
  73.  
  74. /*
  75.  * do the Inverse Fast Fourier Transform
  76.  */
  77. void FFTwrapper::freqs2smps(FFTFREQS freqs,REALTYPE *smps){
  78.     tmpfftdata2[fftsize/2]=0.0;
  79. #ifdef FFTW_VERSION_2
  80.     for (int i=0;i<fftsize/2;i++) {
  81.     tmpfftdata1[i]=freqs.c[i];
  82.     if (i!=0) tmpfftdata1[fftsize-i]=freqs.s[i];
  83.     };
  84.     rfftw_one(planfftw_inv,tmpfftdata1,tmpfftdata2);
  85.     for (int i=0;i<fftsize;i++) smps[i]=tmpfftdata2[i];
  86. #else
  87.     for (int i=0;i<fftsize/2;i++) {
  88.     tmpfftdata2[i]=freqs.c[i];
  89.     if (i!=0) tmpfftdata2[fftsize-i]=freqs.s[i];
  90.     };
  91.     fftw_execute(planfftw_inv);
  92.     for (int i=0;i<fftsize;i++) smps[i]=tmpfftdata2[i];
  93. #endif
  94.  
  95. };
  96.  
  97.  
  98.  
  99.  
  100.